Effects of Strong Correlations and Disorder in d-Wave Superconductors 
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We use exact diagonalization techniques to study the interplay between strong correlations, su- 
perconductivity, and disorder in a model system. We study an extension of the t-J model by adding 
an infinite-range ci-wave superconductivity inducing term and disorder. Our work shows that in the 
clean case the magnitude of the order parameter is surprisingly small for low-hole filling, thus im- 
plying that mean-field theories might be least accurate in that important regime. We demonstrate 
that substantial disorder is required to destroy a d-wave superconducting state for low-hole doping. 
We provide the first bias free numerical results for the local density of states of a strongly correlated 
d-wave superconducting model, relevant for STM measurements at various fillings and disorders. 

PACS numbers: 75.10. Jm,05.50.+q,05.70.-a 



The combination of strong correlations and reduced 
dimensionality makes the theoretical understanding of 
high-T c superconductivity very difficult and consensus 
on its origin has not been reachedjii 2 -'^ Very recently, ex- 
periments using local probes, such as scanning tunneling 
spectroscopy (STS) and scanning tunneling microscopy 
(STM), have shown that the doped cuprates are highly 
inhomogeneous^ (for a recent review, see Ref.H). Specific 
aspects of high-T c superconductivity, such as the robust- 
ness of the tunneling spectrum with respect to disorder,— 
emphasize its contrast with more conventional disorder 
sensitive^. BCS-type superconductivity. 

In this work, we probe the effects of strong correlations 
and disorder in d-wave superconductors derived from a 
Mott insulator. We introduce and study a generalized 
model derived from the t-J model 8 - in which a supercon- 
ducting (SC) ground state is argued to be inevitable. We 
consider a Hamiltonian TL = Ht.j + H-d + T~L ra ndom, with 
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where the sum runs over nearest-neighbor sites, a = 
(f, I), and standard definitions for the projected creation 
and annihilation operators are employed^ 

The t-J model is microscopically justified from cither 
the one-band Hubbard model by a large U expansion or 
more generally from a reduction of the three-band cop- 
per oxide model to a single-band model^ with a greater 
freedom for the parameter ratio J/\t\. While a mean-field 
theory (MFT) for the doped t-J mode l 10 ' 11 gives a d-wave 
SC ground state, it is an uncontrolled approximation. In 
the MFT, there are often states with other broken sym- 
metries in the proximity of the superconductor that can 
be missed. Most importantly, in view of the very strong 
constraint of single occupancy in the model, we expect 
significant quantum fluctuations, and the MFT cannot 
handle these precisely. Therefore it is not clear that the 
t-J model does have a SC ground state for the ranges of 
parameters studied. In order to precipitate a SC starting 



state, we add an attractive term 
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where A = (Ai,i+st 
L is the number of lattice sites. This is an infinite-range 
term of the type that BCS considered in their reduced 
Hamiltoniaiifi 2 - while building in the d-wave symmetry 
of SC order. We have also considered imposing an ex- 
tended s-wave symmetry, where the results are quali- 
tatively quite different and will be reported elsewhere. 
Within MFT, this model leads to the same d-wave state 
as found from the t-J modeh 10 ' 11 Our model is presum- 
ably a superconductor for any Xd ~ O(l) in the thermo- 
dynamic limit, and for sufficiently large Xd, for any rea- 
sonable finite cluster. Notice that we have sidestepped 
the issue of the "mechanism" of superconductivity, which 
cannot be settled with studies of the kind undertaken 
here and focus instead on the nature of the state so pro- 
duced. We argue below that despite the infinite-ranged 
nature of Hd, strong correlations produce a non- mean- 
field-like state; this state has an unexpectedly small order 
parameter (OP). 

Finally, we consider a quenched random disorder term 
of the form 
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where the Ei 's are taken randomly from a uniform distri- 
bution between [— r, T]. The full Hamiltonian Eqs. ((T|)- 
(|3"J) thus describes an inhomogeneous strongly correlated 
superconductor. In our study, we use numerical diago- 
nalization of clusters with 18 and 20 sites. The dimension 
of the largest Hilbert space diagonalized here is ~ 10 8 . 

In our model, we are interested in understanding how 
the evolution into the SC state occurs, as A^ is turned 
on. Towards this end, we show in Fig. [T] the derivative 
of the energy (main panels) and the energy itself (in- 
sets) as Xd is increased, for different fillings of the 18 and 
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FIG. 1: (color online). Energy (insets) and its derivative 
(main panels) vs \d. (a) Low density of particles (less corre- 
lated) and (b) low density of holes (strongly correlated). N 
is the number of electrons, and J = 0.3. All energies are in 
units of t. 



20 site clusters. For low fillings of electrons [Fig. [Ha)], 
we find that for certain cases of the number of particles, 
level crossings occur, as signaled by a jump in the en- 
ergy derivative, indicating a change of the symmetry of 
the ground stated A similar jump is seen for all fillings 
between N = 4 and N = 10 in the 18-site cluster [not 
shown in Fig. [Ha)] . On the other hand, we find that for 
low-hole fillings in the 18 and 20 site clusters, the energy 
derivative is continuous with A^, suggesting a particu- 
lar compatibility between the <i-wave order and the t-J 
model. It is interesting that there is evidence of this com- 
patibility from high-temperature expansions^ and exact 
diagonalization studies^ which are also unbiased such 
as the present one. 

One of our diagnostic tools for studying the nature 
of the SC state is the d-w&ve pair density matrix Py = 
{D\Dj)^ Its largest (Ai) and next largest (A2) eigenval- 
ues are computed and their ratio R(> 1) is monitored. 
This ratio is an effective probe of the order, both for 
clean and disordered superconductors. Taking the ratio 
eliminates uninteresting normalization effects related to 
the change in the particle density, etc. This procedure, 
for example, eliminates the expected diminishing of all 
the eigenvalues of Py as the hole doping decreases (due 
to Gutzwiller correlations). It is thus constructed as a 
pure number. For a SC ground state, it is expected to 
scale for large L like R ~ ^f 2 L + $, where VP is the di- 
mensionless O(l) OP (Ref. [H) and $ - 0(1) represents 
the depletion (i.e., spillover) from the condensate. This 



FIG. 2: (color online). Ratio R — A1/A2 with increasing 
Ad. (a) Low density of particles (less correlated) and (b) low 
density of holes (strongly correlated). In the insets, we com- 
pare the infinite-range model Eq. ] with a short-range ver- 
sion of Eq. ([2} [taking i = j and dropping the L in the de- 
nominator] up to larger values of A^. N is the number of 
electrons, and J = 0.3. 



depletion occurs due to repulsive interactions, i.e., strong 
correlations. In the parallel case of a Bose system with 
N b bosons^ we expect Ai ~ 0(N b ) while A 2 - 0(1). 

In Fig. [51 we show R = A1/A2 for different fillings of 
the 18 and 20 lattices as a function of Ad- By comparing 
the low electron filling [Fig. [2ja)] to the low hole doping 
case [Fig. |2£b)], one can gauge the effects of correlations 
for overdoped (a), and optimal or underdoped cuprates 
(b). For the lowest electron densities (N = 4), R reaches 
very large values (R > 30) and decreases as the density 
is increased, up to around R ~ 6 for ten electrons. We 
notice that sometimes for low electron filling one needs 
Arf to exceed a critical value before R starts increasing, 
e.g., N = 6 in Fig.[^a). This is a signature of a quantum 
phase transition into a SC state, and occurs in many but 
not all instances. The transition point coincides with the 
jump seen in the derivative of the energy in Fig. [1] Taken 
together, these confirm that the new ground state has a 
different symmetry than the ground state of the plain 
vanilla t-J model. 

On the opposite end, for low-hole doping [two and four 
holes in 18 and 20 sites in Fig.^b)], one can see that R 
increases continuously with Ad, i.e., no abrupt transition 
occurs. Figure [2{b) also shows that in that regime R in- 
creases very slowly with A^ and does not exceed R = 2 for 
Ad < 1. We can interpret this as a small value of 'J? and 
a large value of $ as defined above, implying a large de- 
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FIG. 3: (color online). Ratio R — A1/A2 with increasing 
disorder (r). (a) Low density of particles (less correlated) 
and (b) low density of holes (strongly correlated). In the 
insets, we show the separate values of Ai (upper plots) and 
A2 (lower plots), corresponding to the results shown in the 
main panels. Here J = 0.3 and Ad = 1. These results were 
obtained averaging over ten different disorder realizations. 



pletion of the condensate. This shows, that even in our 
infinite-range model, the SC OP at low doping is very 
strongly depleted. However by no means should we un- 
derstand that superconductivity is weaker in that regime. 
From studies of purely bosonic systems, it is known that 
due to strong correlations, the superfluid (SC) fraction 
(i.e., p s ) can be much larger than the condensate frac- 
tion (i.e., the OP "J) (see, e.g., Ref. [I||). Correlations 
also tend to make superfluidity (superconductivity) more 
stable against perturbations. 

We now compare our results for the infinite-range 
model of Eq. ||3J) to those produced by the more stan- 
dard short-range case [Eq. ([2]) , for i = j and no normal- 
ization by L in the denominator] used in the literature 
dealing with the Hubbard model. The insets in Fig. [2] 
show that while in the infinite-range model R saturates 
with increasing A^, in the short-range model R attains a 
maximum value for Xd ~ 1 and then decreases towards 
unity. The latter occurs because for A^ ;§> 1 one produces 
localized pairs, i.e., there is no long-range coherence. 

Turning to disorder, in the main panel in Fig. [3J we 
show how R evolves with increasing disorder for A c ; = 1. 
(Those results were obtained averaging over ten different 
disorder realizations.) Here one can see that disorder 
produces a very large reduction of Ai / A2 for low electron 
filling, i.e., disorder has a very large impact on the SC OP. 
For the case of low-hole density, the relative reduction of 



A1/A2 is much smaller, i.e., as Xd had a small effect in 
increasing R, so is V having a smaller effect in reducing 
it. 

From the main panels in Fig. [3J we see that the effect 
of disorder in the SC state is always to make R decrease. 
It is of considerable interests to understand what hap- 
pens to Ai and A2 separately as T is increased. Results 
for these quantities are presented in the insets in Fig. 
[3j There one can see that Ai behaves qualitatively very 
differently between low-electron fillings and low-hole fill- 
ings. In the first case Ai exhibits a very large reduction, 
which points towards the destruction of superconductiv- 
ity On the other hand, for low-hole doping, Ai is almost 
unaffected by the increase of disorder and can even be 
enhanced, as shown for 14 particles in 18 sites. Unex- 
pectedly, the reduction of R in this case is related to an 
increase of Aj. This increases points towards a slower 
decay of Pij when disorder is increased. This suggests 
the possibility of an emerging algebraic long-range order, 
producing a different signature in the density matrix than 
the case of standard LRO. For example, in the 2D XY 
model below T c , or in the ID Hciscnberg antifcrromag- 
netic ground state, there is no true LRO, but several of 
the largest density matrix eigenvalues scale as L v , with 
77 < 1. The system sizes we treat here are too small 
to make definitive statements. However, it is interesting 
to note that for low-hole doping the behavior is quali- 
tatively different from the low electron filling, in which 
the largest eigenvalue exhibits a large decrease with in- 
creasing disorder. An analysis of the data for the s-wave 
superconductor studied in Ref. Qj^ exhibits exactly the 
latter behavior, in contrast to the one we see for the SC 
t-J model in the low-hole doping regime. Our results 
therefore suggest unusual power-law type superconduc- 
tivity in the presence of disorder close to half filling. 

In order to make connection with experimentally mea- 
surable STM curves, we show in Fig. [4] the local density 
of states of the L = 20 site cluster for two different fillings 
in the presence of disorder and A^- Figures HJa) andfjjd) 
correspond to fillings where the ground state of the plain 
t-J model (Xd = 0) is adiabatically connected to the SC 
ground state at finite A^. In the presence of disorder, the 
density of states is similar to the one reported previously 
for the translational invariant t-J model with L = 16 
sites^ These curves display a striking asymmetry be- 
tween adding a particle and taking out a particle, and 
the evolution of this asymmetry with doping is similar to 
that of the clean t-J model. 

Adding the SC term (A^ > 0) to the disordered sys- 
tem opens a gap. This can be clearly seen in Figs. H{b) 
andfjje). Our system sizes are too small to see the V 
shape expected for a d-wave superconductor, i.e., we see 
a real gap. As disorder is increased, Figs. (He) and|4jf) 
show the reduction of the gap. From the results shown in 
Figs.Hl we see that the SC gap closes only for a substan- 
tial disorder (r > 2Xd). Our calculations therefore also 
sheds light on this aspect of the STM spectra, namely, 
the robustness against disorder. 
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FIG. 4: (color online). Averaged density of states [JV(w)] 
for two different fillings (N = 8 and N = 16) of the 20 site 
cluster. We show results for: T = 1, Ad = (a) and (d); 
r = 1, Ad = 1 (b) and (e); and r = 2, X d = 1 (c) and (f). 
N(ui) was computed as the average over different lattice sites 
and over two different disorder realizations. 



In conclusion, we have presented and studied a vari- 
ant of the t-J model, with an infinite-range d-wave su- 
perconducting term. We have shown how the energy, its 
derivative, and the d-wave superconducting order param- 
eter evolve with increasing the strength of the supercon- 
ducting term. In addition to discontinuities in all the 
above quantities for low electron densities, we find a se- 
vere reduction of the magnitude of the order parameter 
at low-hole filling. This is a signature of strong quan- 
tum fluctuations near the Mott insulator. In relation 
to current STM experiments, we find that superconduc- 
tivity survives considerable disorder close to half filling. 
The local density-of-states curves yield bias free (i.e., non 
variational) results for a strongly correlated d-wave su- 
perconductor in the presence of disorder and provide a 
picture of the large energy scale structure of this impor- 
tant object. 
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